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1. Introduction 

The use of holographic methods to explore strong coupling in gauge theories has 
yielded a particularly fruitful interaction between string theory and low temperature 
physics. The typical set-up uses the concept of gauge/gravity duality [Q in which a 
classical gravitational system with negative spacetime curvature has, on its boundary, 
equivalent degrees of freedom to a strongly coupled gauge theory. Temperature is 
gravitationally introduced into these systems by adding a black hole in the bulk 
spacetime, and different holographic dual theories can be constructed by having 
additional bulk fields (see [0 for reviews of this approach). 
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The standard application of gauge/gravity duality is the adS/CFT correspon- 
dence, which yields a generally scale invariant boundary theory, however, more re- 
cently, attention has focussed on systems having more general scaling properties, 
such as non-relativistic field theories, 0, or, pertinent to this investigation, a gen- 
eral dynamical Lifshitz scaling, z: 

t^\H , x' ^ \x' , r^ r/\. (1.1) 

In order to produce such a dynamical scaling, the spacetime metric must be posited 
to have the following form 

ds^ = L^ (r^'de - ^ - r''dx4x'\ , (1.2) 

which explicitly respects the scaling ( |1 . 1[ ) . In this metric, not only the asymptotics, 
but the full spacetime has the required scaling symmetry. Clearly, such a spacetime 
requires a matter content to produce this asymmetry, and this was first set out in 
the paper of Kachru et al. [Q , in which charges and fiuxes of topologically coupled 
gauge fields provided the necessary scaling. This theory is in fact on-shell equivalent 
to a somewhat simpler massive vector theory |^, although the r — )■ singularity of 
these spacetimes exhibits certain pathologies ^. 

As with any holographic theory, although we can explore empirical simple mod- 
els, in order to have confidence that there is indeed a holographically dual field theory 
we should be able to construct a qualitatively similar "top down" theory with Lifshitz 
scaling within string theory. After initial halted progress, string theory embeddings 
of Lifshitz geometries with dynamical exponent z = 2 were found in J^, |[, by making 
a consistent massive truncation of type IIB supergravity to a lower dimensional the- 
ory resembling the phenomenological construction of [Q. Soon after, a method for 
constructing Lifshitz spacetimes within string theory for arbitrary scaling exponent 
z > 1 was put forward in P|. In this approach, the Lifshitz space is constructed 
from a simple fiux compactification of Romans' gauged supergravity in five and six 
dimensions, [|10], [ll[], generalising the classic adS compactifications of those theories. 



The lower dimensional supergravity theories can be obtained by dimensionally re- 



ducing type IIA or IIB supergravity, as shown in [12, 13 1, and any solutions can 



immediately be uplifted to ten dimensions^. Further Lifshitz and AdS solutions in 



gauged supergravity and string theory have been also studied in |T5 . 

In order to explore physical dualities, we need to be able to set our system at 
finite temperature, in other words, we need to introduce a black hole to our space- 
time. Black holes in asymptotically Lifshitz spacetimes were initially hard to build. 
However by now several such solutions have been found in simple phenomenological 
models, starting with the numerical work of [TH] . By engineering a matter or gravity 



^See mM for other examples of non-relativistic solutions in massive type IIA supergravity. 



-2- 



content to source the desired geometry, some analytical solutions have also been con- 
structed. Overall, Lifshitz black hole solutions in phenomenological models include 
numerical and analytic studies, fixed as well as arbitrary critical exponents, hori- 
zons with various topologies, extensions to other dimensions, higher-order theories 
of gravity, and Brans-Dicke models [|17|, 0, |1^, |2^ . 



As to embedding Lifshitz black holes into string theory, this can now be done fol- 
lowing on from the string constructions of pure Lifshitz geometries discussed above. 
Recently, numerical string Lifshitz black holes with dynamical exponent z = 2 were 
presented in I^Tj. As in 0, ^, their method was to identify a consistent massive trun- 



cation from type JIB supergravity to a lower dimensional model resembling previous 
phenomenological constructions. In the present paper, we construct string Lifshitz 
black holes with general dynamical exponent z > 1, generalising the lower dimen- 
sional supergravity/type IIA Lifshitz solutions, which were found in by deforming 
adS solutions. Alongside the asymptotically Lifshitz black holes, we study related 
asymptotically adS black holes; thus we are able to draw on the intuition gleaned 
from the latter as well as identify which properties belong uniquely to the Lifshitz 
case. 

With holographic condensed matter applications in mind, our interest is in planar 
black hole geometries, whose boundary field theory propagates in flat 2+1 spacetime, 
and moreover we consider static geometries corresponding to equilibrium phases. 
Naturally, the black hole solutions are not so simple as their pure adS or Lifshitz 
cousins, with the exception of adS-Schwarzschild. By exciting the supergravity fields 
about this latter background in the probe limit, we learn about the charges inherent 
to our system. Further progress can be made by expressing the supergravity field 
equations as an autonomous dynamical system, whose fixed points are pure adS or 
Lifshitz. By perturbing close to the fixed points we can understand in detail how 
general interiors, including black holes, can fiow to the adS/Lifshitz asymptotics. 
Indeed, this method allows us to analytically characterize all the possible asymptotic 
behaviours for static adS and Lifshitz black holes for our theory and moreover helps 
to numerically integrate to the full black hole solutions. 

As might be expected, since our string/supergravity setup contains more degrees 
of freedom than the simple phenomenological models, the black holes have a rich 
structure. In the end, the black holes we find necessarily have some non-trivial scalar 
field, and aside from the horizon size, are characterized by two parameters, which 
can be interpreted as a form field charge and scalar charge. Thus we can begin 
to explore their configuration space, how the field profiles and thermodynamical 
properties change as the charges and dynamical exponent z vary. 

The paper is organised as follows. In the next section, we present the six di- 
mensional supergravity theory where our black holes will be found, its pure adS and 
Lifshitz solutions and the general planar geometries that we will study. In section 3 
we proceed to analyse this setup in detail, starting with some approximate analytic 
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solutions, both through the whole spacetime but close to the adS-Schwarzchild black 
hole in the probe limit, and in the far field limit close to the asymptotic adS and 
Lifshitz geometry. We then build upon these results in section 4, to find numerical 
solutions describing adS and Lifshitz black holes, uplift them to type IIA super- 
gravity, and study their behaviour. In section 5 we briefiy discuss thermodynamical 
properties of the black holes and we conclude in section 6. We give the details of the 
dynamical system we use to solve the supergravity field equations in appendix A. 
Finally, as an aside, in appendix B we identify some exact analytical black hole solu- 
tions in a generic dilatonic model, which could plausibly be related to a supergravity 
theory. 



2. The System 

In this section, we introduce the supergravity theory that will be the subject of the 
paper, present its pure adS and Lifshitz solutions, and propose the general Ansatz 
which we use to find planar black holes that asymptote adS and Lifshitz geometries. 
We consider six dimensional A/" = 4 gauged supergravity, first presented by Romans 
in 



I0| . This theory can be obtained from a consistent truncation of massive Type IIA 



supergravity, and thus the solutions of the six dimensional theory can be uplifted to 



solutions in string theory [|12| . In |Q it was found that the field content and couplings 
of this theory admit Lifshitz solutions. 

The bosonic field content of 6D Romans' supergravity consists of the metric, 
qab, a dilaton, 0, an anti-symmetric two- form gauge field, Bab, and a set of gauge 
vectors, {A)^ ,Aa) for the gauge group SU{2) x f/(l). The bosonic part of the action 
for this theory is 






— - £^ ^^ Bab ( J^cdJ^ef + ^^BcdJ^ef + -;^BcdBef + Fc'd-^ef 



1 ABCDEF n .^ I T-^^T-__ I _D„_Tr__ i ^ D„_R__ , l?» fW 



f2.11 



where g is the gauge coupling, m is the mass of the 2- form field Bab, J^ab is a 
U(l) gauge field strength, F^*^ a nonabelian SU(2) gauge field strength, and T-Lab = 
J^AB + nii?AB- Spacetime indices A, B, ... run from to 5, and e is the Levi-Civita 
tensor density. Notice the presence of Chern-Simons terms in the previous action, 
identified in M as an important ingredient for the existence Lifshitz configurations. 
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Varying the action gives the equations of motion: 

Rab = 2dA<pdB(p + lgABV{(P) + e^^'^ (g^^Gbcd - IqabG' 

°^ = l^+ l]ll''^"'^" + l]l'r^'' (^' + (^'^)') (2.3) 

Vb (e'^^^n^^) = I e^^^^^^UBcGuEF (2.4) 

Vb (e-^-^FW^^) = i e^^^^^^Fi;),Gz.^^ (2.5) 

Vc (e^v^^G^^^) = -me-^^^U^^ - i e^^c:'Di.F [ucdUef + i^<?ii^S) ,(2.6) 
where we have defined the scalar potential function: 

V{ct)) = - [g^e-^^ + 4m(?e-^<^ - m^e-^^^"^) . (2.7) 

Analogously to the Romans' solution adS4 x H2 fl^ , it was shown in that one 
can have a Li4 x H2 dimensional reduction of this 6D supergravity to a 4D Lifshitz 
space with an internal hyperbolic manifold threaded by non-abelian magnetic fiux. 
The solution is given by 

ds^ = L^ (r^'de - r^dxj - r^dxj - ^] - -Adyl + dyl), (2.8) 

V r^ J yi 

where a is a constant, the radius of curvature of the hyperboloid which can be taken 
to be compact (see |2^ for details). The dilaton is also chosen to be constant, = 0o; 



and the field configurations are 

% (2.9) 

The relations between the various constants are somewhat simplified by performing 
the following rescalings 

r- (2-10) 

g = Lge'l"-''^ d = a/L m = Lm g-^-^"/^ . 



Equations (p^ ) to ( p.6|) then reduce to a simple set of algebraic equations with the 
following general solutions 



?2 , ^2 o ^. N ^^ 6 + ^^2^2(2 + 4) 



b' = z-1 g' = 2z{A + z) — 



c2 



Q 



{2 + z){z-3)±2y/2{z + 4) 1 



2 ^ (2.11) 



2z a 
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= Q + SzT2^/2{z + A). 



These define two families of Lifsliitz spacetimes, one for each branch of the square 
root in (|2.11|) . The requirement that b is real restricts z to z > 1 and for the lower 
sign choice, for q to be real one finds that z must be greater than approximately 
4.29. 

In addition to these Lifshitz solutions, the system also allows for an independent 
one parameter family of adS solutions when z = 1 and 6 = 0. These solutions can 
be parametrized by either g or rh, the latter case giving 

m2 + 6 ., 5m^ - 36m^ + 36 1 _ 5m* ~ Vlrn^ + 36 



9 



q 



im^ 



(2.12) 



2m ' "" IQm'^ 

In this case the requirement that q is real implies that m G [a/|, V^]- When m = 

a/IO — 2 the adS and (upper sign) Lifshitz solutions touch at ^ = 1. Figure |I| 
represents the values of z and m that these Lifshitz and adS solutions can take. 




Figure 1: Plot showing the values of z and rh that the Lifshitz and adS solutions can 
take. The horizontal dashed line indicates the adS solutions, with z = 1. The black line 
corresponds to the Lifshitz solutions with the upper sign choice in ( |2.11] ) and the grey to 
the lower sign choice. Notice that adS and Lifshitz solutions meet at m = yTo — 2. 



In what follows, we continue to analyze adS and Lifshitz solutions, generalizing 
the previous discussion to the case of asymptotically adS and Lifshitz space-times. 
Insights acquired when analyzing asymptotically adS configurations will be of great 
help when considering asymptotically Lifshitz space-times. 

2.1 General planar spacetimes 

Our main aim is to characterize asymptotically adS and Lifshitz black hole solutions. 
We look for solutions which respect the planar symmetry and static nature of the 
metric (|2.8| ), meaning the alterations need only have radial dependence, = 0(r) 
and 

ds'' = L^ [e^^^dt^ - e^^^^dx" - e'^^^'Ur''] - e^^^'-^dHl . (2.13) 
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We choose our field strength Ansatze to be 



Q 



Fr = L^O(r) F(3) = -L b = —e 

1/2 ^ 



-V2<t)o 



P(r) 



which gives the gauge equations: 



^2c+2h-f-d-V2(l> 



Q 



qe 



-V2<t>0 JD' 



(2.14) 



(2.15) 



rf-d+2h-2c+2V2c^^-V2<Pop'y ^ ^2j2^-V24>op^f+d-2c+2h-V2^ ^ 4^2^g ^ (2.16) 

Integrating (|2.15|) and noting that Q — > as P — ?■ from ( |2.16|) , we obtain: 
hence there is a single equation of motion for the gauge fields: 

ff~d+2h-2c+2V24>pi\ ^ j2p^f+d-2c+2h~V24> f^2 _^ ^2^2V24>-ih\ ^2 Ig) 

The remaining equations are: 

n/2 f \/9V_\ _ p2 -2v^«!.o-V20-4c ( 2 2V2<t>-ih _ ^A , -^ 2V2(<}!.-0o)-4c-2<J 

_^2g-v^,^-4^ - - (^g'e^^ - Amge-^'^ + Sm^g-^^^-^) (2.19) 

2 ( \/9J_\ _ p2^-2v^9io-v^0-4c / q„2 2v^?i-4h , "^ 1 , -^' ^2V2(0-9io)-4c-2d 

^^2g-y2</,-4h ^ 1 j^^2gV20 ^ 4m(7e-^'^ - m^e-''^^^ (2.20) 



2g-2v^</.O-V20-4c I _ 2^2V2<j)-4h 



±fV9±\ = P 

^^2^~V2<p-Ah ^ 1 (^g2^V2<^ ^ 4mge-^'f' - n 

_2_ [ VdJ^] _ p2-2V24>o-V2^-4:C ( _^2„2V2^-4h . ^\ . 



?/2 



^ „2V2(,^ 

2^2^^ 


-•/-o) 


-4c-2d 


^2g-3v/2<^A 
/ 






p/2 


0o)- 


4c- 2d 



(2.21) 



2L2 



-Sg^e-v^^-^'^ + - [g\^^^ + 4m(7e-^<^ - m^e-^v^-^) - 26-^'^ , (2.22) 



together with the first integral of the Einstein equations: 

p/2 



2/'c' + 2fh' + Adh' + c'^ + /i'^ - <P" - — e 



2V2(</.-0o)-4c 



(2.23) 



2„2d 



L'e 



_g-2/i _ ^2^-V2<P-4h _ p2^-2V2<Po~V2<t>-'ic I ^2^2V2^-4h _^ ^^ j _)_ y (^^\ 



In analysing the solutions of these equations, it is particularly useful to consider 
the equations of motion from a dynamical systems perspective. The exact Lifshitz or 
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adS geometries correspond to fixed points of tfie dynamical system, and a perturba- 
tive analysis around tlie fixed points indicates the flows of general interior solutions 
to the asymptotic Lifshitz or adS geometry. This general methodology was used in 
[p3[] to explore flows between exact Lifshitz and adS solutions. 

Although the system ( |2.18| - pr2^ appears to be nine dimensional, there is a re- 
dundant gauge degree of freedom corresponding to a rescaling of the coordinates, 
and also the Bianchi identity ( ^.231 ), which reduces the order of the system to seven. 
Perturbing around a fixed point solution therefore will give a seven dimensional solu- 
tion space (some of which will correspond to unphysical singular solutions), spanned 
by the eigenvectors of the perturbation operator around the critical point, with a 
radial fall-off given by the corresponding eigenvalue A. In appendix ^ we provide 
the details of a dynamical systems analysis of these equations of motion, and the 
derivation of the eigenvalues. The eigenvalues will provide crucial information about 
the physical charges that characterise a given solution. 



3. Analytic results 



In order to develop a general understanding that will be later used to determine 
numerical black holes configurations, it is useful to analytically explore the allowed 
asymptotic behaviour of adS and Lifshitz solutions for our theory. In this section, 
by implementing the dynamical system analysis developed in appendix 0, we will 
first identify perturbative solutions which describe flows towards an adS boundary 
at large r. Among these is the exact adS-Schwarzschild solution, and we analyse lin- 
earized solutions about this background, throughout the space-time from horizon to 
boundary. We then study flows towards a Lifshitz boundary at large r, for arbitrary 
dynamical exponent z. In this way, we are also able to observe how asymptotically 
Lifshitz geometries for 2; > 1 reduce to asymptotically adS geometries at 2 = 1. 

For convenience, we choose the radial coordinate to correspond to the area gauge 
as in (|1.2|) , i.e. cir^ = logr, and rewrite our metric and scalar functions in terms of 
deviations from the known Lifshitz background: 



{r) = (f)o + (p{r)/v2 , /i(r) = loga + |lniJ(r) , P{r) = r'^p{r) 

(3.1j 



/(r) = 2; Inr -|- ilnF(r) , d{r) = — Inr — | In D{r 



This gives us the field equations 
{r'+^VFD Hip')' q^e-^ 



r.2 + 1 



H^/FjD 



{r'+^VFDH')' 
r'+^H^F/D 






2 / g^e"^ m^e"^ 



p 



rp 



dV 



2De'-(p+^] -— (3.2) 



-P 



2 , q^e'P rn^e-'P 






H^ 



4 



+ 2De'^'P I p + 



rp 



^2+3 



..z+l 



^D/FHF')' 
H^F/D 









(3.3) 



H"^ 



-2{l + z)De''P \p + 



rp 



{2r^+^HVFD)' fe' 



P 



q^e^ 3m2e-^ 



r^+^H^/F/D H^ ' \ H^ 

together with the gauge equation 



+ 



+ {l-z)V{^) 
2De'^{p+^4- 



(3.4) 
V{ip), (3.5) 



VfD e^'P H {2p + rp') 



r 



-^pHe-^^/F/D(m' + ^^ 



(3.6) 



and the first integral 
rF' 



D 



2z + l + 



,rH' r'^H'F' r'^H''^ 

iz + 2)^^ + ——— + 



^2^/2 



H 

Vicp) 



2HF 



1 q^e-^ 



am 



H 



AH'^ 2 

P 



P+"-{\ e^- 



2 / q^e^P rfi^e-'P 



H^ 



where 



^(V^) = ^ (^^e'^ + Amge 



"f _^2^-3</p\ 



m e 



(3.7) 
(3.8) 



3.1 Anti-de Sitter solutions 



We start with a detailed analysis of the adS branch of solutions to our system. These 
configurations are easier to study than Lifshitz space-times, yet they give insights 
into what charges are inherent in the system, and what aspects of black hole solutions 
are uniquely Lifshitz. In the adS case, of course, we already know an analytic black 
hole, the adS-Schwarzschild solution: 



H 



ip = p = {). 



r 



(3.9) 



From the point of view of the dynamical system discussed in appendix ^ this cor- 
responds to the nonlinear evolution into the interior of the asymptotic eigenvalue^ 



^Asymptotic eigenvalues A^ control the fall-off r '^' of a given function at large 
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A = —3 (the 1/r^ fall-off), which has an eigenvector with components only in the 
directions corresponding to deformations of the 4D geometry, as in eq. ( p.9|) . 

It is useful however to continue with a more general analysis of perturbative 
solutions which asymptote adS, as this will enable a more direct comparison with the 
Lifshitz case. A general analysis of the perturbations around the adS fixed point (see 
appendix ^ yields the eigenvalues plotted in figure Q. As explained in detail in the 
figure caption, these eigenvalues give the exponents of r in the asymptotic solutions 
for the various fields, as function of rh. Each field has a fall-off of a pair of exponents 
which are symmetric about —3/2, and whose coefficients can be interpreted as a 
source and operator in the boundary field theory. The pure 'black hole' mode, 
A = —3 is in this sense dual to the zero mode which takes us along the parameter 
space of adS solutions. 



Re{A,) 




Figure 2: Plot showing the real part of the field exponents as they asymptotically approach 
an adS space-time, as a function of the parameter m. Each pair of exponents sums to —3. 
The solid black lines correspond to switching on a B-charge in the interior spacetime, 
and the solid horizontal grey line is the black hole solution. A combination of the dotted 
and dashed black and grey lines corresponds to switching on the dilaton and breather 
modes. (The joining of the dotted black lines for small m indicates that the exponents 
turn complex.) Notice that at m = ylO — 2, which is the value of m at which the adS fixed 
point solution is equivalent to a 2 = 1 upper branch Lifshitz solution, there is a degeneracy 
in the eigenvalues: the mass and charge deformations have the same fall-off. 



A combination of the dotted and dashed black and grey lines in figure ^ cor- 
responds to switching on the dilaton and breather modes, keeping P equal to zero. 
The solid black lines in the figure are associated with the turning on of the 2-form 
charge P only. At m = a/IO — 2, the value of rh for which the pure adS and Lifshitz 
solutions coincide (since z = 1), and the charge and mass deformations of adS be- 
come degenerate. Since the P-equation decouples at leading order, it is not difficult 
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to extract the charged perturbations of adS: 

H = F = D = 1, ip = 0, p = po-^. (3.10) 

The po deformation corresponds to the zero mode (the upper sohd black hne at 
m = a/IO — 2 in figure 0) which moves the solution onto the Lifshitz branch whereas 
the ps deformation identifies the pure charge eigenvector. 

As we learn in what follows, perturbations of asymptotically Lifshitz configura- 
tions have different eigenvectors and eigenvalues, conveniently parameterized by the 
quantity z. On the other hand, when m = vlO — 2 (corresponding to 2; = 1) the 
Lifshitz deformations and adS deformations should coincide, as there is a degeneracy 
of the eigensystem at this point: the corresponding Lifshitz deformations are a com- 
bination of the pure mass and pure charge asymptotically adS solutions. In Section 



3^ we analyze the corresponding Lifshitz asymptotic solutions and how they join 
the adS ones at 2; = 1. 

Turning now to black hole solutions, we next explore linearized solutions for 
scalar and gauge charges around the known adS-Schwarzschild black hole background 
of eq. (13.91). We stress that the eigenvalue analysis discussed above refers to pertur- 
bations of the full system around the adS background, and as such describe how the 
geometry and fields asymptote the adS boundary at large r. In contrast, in the next 
two subsections we seek linearized solutions of either the 5— field, or the dilaton and 
breather mode, around the black hole background. These, at leading order, do not 
include perturbations of the black hole geometry, but they are solutions of the scalar 
or gauge fields for the full range of the space-time from the horizon to the boundary. 
Our linearized solutions, therefore, should asymptote one or more of the eigenvalue 
solutions near infinity, at least for the fields we are perturbing. On the other hand, 
they are more informative since they provide some further understanding on how 
field configurations behave in the black hole background. Our aim here is to identify 
the eigensolutions of the full system, so as to have a physical interpretation of the 
various eigenvectors which we can then use to understand the Lifshitz system. 

3.1.1 i?— charge 

We start by allowing p to vary, but keeping (p and H constant over the black hole 
background of eq. ( p.9|) . Since P is related to the electric field and the flux of the 
Bab field of the system we shall refer to these solutions as charged black holes. Note 
that, to leading order, the p— equation (|3^ ) decouples from the other equations and 
around a black hole background is 

^{^+^f)p- C"'^-^-"' )p. (3.11) 

Writing x = {r+/r)^ and definining 

u± = (±V(36 - m2)(m2 - 1) - m)/6m, (3.12) 
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1 - ^ ) (rV)' 



we obtain a closed solution for P = r'^p in terms of hypergeometric functions: 



P{r) = r[2z/_ + 4/3] r[z/+ + 4/3] T[u+] x' 
-T[2i^+ + 4/3] r[z/_ + 4/3] T[uJ x 



2FiK,z/+ + 4/3,2z/+ + 4/3;x(r)] 
"- 2Fi[z/_, u_ + 4/3, 2i._ + 4/3; x{r)] 



(3.13) 



where the constants are chosen to give a nonsingular combination at a; = 1, the posi- 
tion of the horizon. A quick glance at the eigenvalue plot, figure 0, shows that, unless 
m < v^IO — 2, the z/+ hypergeometric function will blow up at infinity. Thus for non- 
singular linearised solutions we take m < vlO — 2 (this does not mean that charged 
black holes do not exist for m > a/IO — 2, simply that they have strong gravitational 
backreaction). Figure ^ (on the left) shows a representative sample charged 5— field 
around an adS black hole with g"^ = 52/5. We will further numerically explore these 
black holes in section ^ 



=^ 0.5 




Figure 3: The left plot shows a perturbation of p for rh = v 10 — 2 whose asymptotic 
behaviour is consistent with ( |3.9D . The right plot shows perturbations of ip (in purple) and 
H (in black) for m = 3/2, which are consistent with { \i.9\) . 



3.1.2 Scalar charge 

From the perspective of the 4D geometry, the breather mode of the hyperbolic geome- 



try appears as a 4D scalar, and indeed the scalar equations, (|3.2| , p.3|) are independent 
of perturbations of the geometry at linear order. Extracting these scalar equations 
gives a second order system: 



C 



6ip 
6H 



■(3m2 - g^)/2 2g2 

2f 2(g2 + 3)_ 



where C is the linear operator 



J.2 ^j. 






dX 
dr 



6{p 
6H 



(3.14) 



(3.15) 
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Diagonalising the matrix on the RHS of ( 3.14|) yields two eigenvalues and eigenvectors 
which correspond to the two pairs of exponents indicated by the grey dashed and 
black dotted lines in figure ^ The fully coupled system also has perturbations of the 
geometry, although the asymptotic exponents indicated in figure |^ represent the fall- 
off at large r. The eigenvalues of ( p.l4| ), however, are linearized solutions around the 
given black hole background for all r. Clearly, examining the large r behaviour from 
figure ^ shows that the dotted grey branch cannot yield a solution which is regular 
at both horizon and infinity: hence these branches have significant backreaction on 
the geometry. However, for m < (6 — ^/6)/V5 the black dotted branch gives the 
regular solution: 



6ip 
6H 



Irh^ ± V36 - 60m2 + 89m4 
5m^ — 6 



X 



r 



(3.16) 



where 



X[x\ = r[2/i_]r[/i+]V+2Fi[;U+,/x+,2/x+;x] -r[2/i+]r[^_]V- 2Fi[/i-,/U_,2/i„;x] 



and 



/i± 



1± v/l + 4A/9 



are given in terms of the eigenvalue 

A = f 3m^ + 36m^ - 36 + (m^ - Q) V36 - QOm^ + 89m'^ 

8m^ V 



(3.17) 
(3.18) 

(3.19) 



As is noted in ||23| the straight part of the dotted black curve, where ^/Q/5 < 
rh < 1.254, indicates an imaginary exponent. This occurs when A < —9/4 in ( p.l9| ) 
and is analogous to a mass violating the Breitenlohner-Freedman bound, p^. This 
fact is indicative of a possible instability, most likely a flow from one adS branch to 
another. Figure ^ (on the right) shows the field profiles for the representative value 
m = 3/2. 

3.2 Lifshitz solutions 

We now apply to Lifshitz configurations the same techniques we developed to char- 
acterize adS solutions. As we will see, the intuition we developed for adS will help 
in characterizing the richer Lifshitz configurations. Unlike the adS case, there is 
no straightforward analytic Lifshitz black hole solution for our system^. That this 
will be the case can be seen by analyzing the Lifshitz fixed point, where all of the 
eigenvectors corresponding to perturbations around the asymptotic Lifshitz geometry 
generically have components in every field. In this instance, there is no pure geo- 
metric deformation to the Lifshitz space, and any deformation necessarily includes a 
scalar and gauge profile. 



^See however appendix [H for examples of exact analytic solutions in similar dilatonic theories. 
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Figure 4: Plots of the real parts of the eigenvalues of perturbations around the Lifshitz 
solution as a function of z. The left and right plots correspond to the upper and lower sign 



choices in ( 2.11 ) respectively. For the upper sign choice, the eigenvalues are all real only 
in the region 5.69 < z < 5.83, whereas for the lower sign choice one finds real eigenvalues 
for z > 16.82. 



The analysis of the Lifshitz point is given in appendix 0, and the system of 
eigenvalues plotted in figure |. Comparing this plot to the adS case, we see that 
the eigenvalues are symmetric around —{z + 2)/2, and in particular, an eigenvalue 
A = — (2; + 2) exists, which is continuous with the pure black hole adS eigenvalue, 
A = — 3. Indeed, plotting the Lifshitz and adS eigenvalues side by side shows how 
the perturbations around the critical points merge as 2; — )■ 1 or m — )■ vTO — 2. Using 
the intuition obtained from analysing the adS solutions, one might expect that the 
mass and charge perturbations continue from the adS side into the Lifshitz side. 
However, things are not so straightforward, and are actually more interesting. While 
the subspace spanned by the two eigenvectors clearly is the same on each solution 
branch, the eigenvector basis need not be: the adS eigenvectors are either pure charge 
or pure geometry, whereas the Lifshitz perturbations include all fields. A careful 
tracking of the perturbations as 2; — )■ 1 indicates that the Lifshitz perturbations can 
be thought as corresponding to a 7r/4 rotation of the adS perturbations. Thus, the 
"—(2; + 2)" eigenvector is actually a charged black hole, most likely an extremal black 
hole given the combination of adS eigenvectors. 



It is worth exploring in more detail the analytic expansions of the functions, 
especially for the special cases identified above. First of all, by analysing the linear 
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perturbations for the eigenvalue A = — (z + 2), we find the solution (for general z) 



6ip 



2MV^^T 



r.2+2 



1- 



2Mv/7^I 



2(4 + z) (-48 -8z + Uz"^ + Az^) 
+ (136 + 40z - mz^ - 21/ - 2/) 

v/2(4 + 2) (96 - 352^ - 228^^ - 22^^ + Az^) 



D = l + 



(2 + z) r^+2 

+ (-272 + 952z + 760/ + 146/ - 11/ - 11/ - 2/) 



.z+2 



^/2{4: + z) (-48 + AOz + 38z^ + 4/) 
(136 - 96z - 124/ - 19/ + 7/ + 2/) 



H 

P 



1 + 



2Mv/^^I 



,2+2 



y^^T + 



2M 

»,2 + 2 



2(4 + 2) (-48 + 8z + 2/) + (136 - 8^ - 12/ + 5z^ + 2/) 
v/2(4 + z) (-48 - 82 + 22/ + 6/) 



+ (136 + AOz - 64/ - 23^^ + 7/ + 2z^) 



(3.20) 



where M is some integration parameter. This solution is valid at first order in pertur- 
bations around the pure Lifshitz solution with arbitrary dynamical exponent z, and 
makes manifest that all fields are normally switched on for Lifshitz configurations. 
We additionally checked that the system of equations can be solved also at next to 
leading order, providing corrections to the above profiles that scale as M^/r^^^"''^): 
the expression for the 2;— dependent coefficients is however too long to be presented 
here. Notice that the limit z — )■ 1 is well behaved, and leads apparently to the 
perturbative pure charge adS solution we discussed in equation (|3.10|) . 

We can also characterize the solution corresponding to the black line in the left 
panel of figure ^, that joins with the grey line when 2; — )■ 1. In order to do this, we 
expand the eigenvalues tending to the A = —3 near z -^ 1 to obtain: 



Ai = -3-(^-l) 

A2 = -3 + -^(260v^ - 701)(2 -l) + OUz- 1)2] 
189 



(3.21) 
(3.22) 



The first is just the eigenvalue A = —z — 2 rewritten as an expansion around z = 1, 
with a corresponding solution, ( p.20|) , that can also be expanded near z = 1. The 
second eigenvalue, A2, is instead that of the black curve in figure H, and we can 
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similarly determine its corresponding solution. To leading order in 2; — 1 we find: 






Pi = V^^+^ , P2 = V^^+^ (3.23) 



for some integration constant /x. Notice that both these eigenvector solutions ap- 
proach the pure charge adS solution for in the limit z -^ 1. On the other hand, for 
small values of {z — 1) they differ (by an identical amount) only in the metric compo- 
nents F and D (up to corrections suppressed by powers of \/z — 1). The difference 
in the eigenvectors is consequently due to contributions of pure geometry to the Lif- 
shitz configuration, that we know corresponds to the adS-Schwarzschild eigenvector 
in the pure adS case. In this sense, we can regard the two eigenvectors ( |3.23D , in the 



limit of small (2; — 1), as if each forming a ir/A degree angle with the eigenvectors of 
asymptotically adS configurations. 



4. Numerical Black Hole Solutions 

Having developed an analytical understanding of the asymptotic properties of the 
solutions of the system, we now present some numerical solutions to the fully coupled 
system of field equations. In this section, we begin by computing and analysing adS 
black hole solutions, followed by their Lifshitz generalizations, and then uplift the 
solutions to type IIA supergravity. 

To obtain black hole solutions for our system we must ensure that our boundary 
conditions are consistent with the nature of the near horizon region of a black hole 
spacetime. These conditions will be the same irrespective of whether we are interested 
in asymptotically Lifshitz or adS black holes. Assuming that the horizon is non- 
degenerate, we wish the gu component of the metric to have a simple zero and the 
Qrr to have a simple pole at r = r_(.. Checking that the matter and metric fields 
and the energy momentum tensor are regular at the horizon imposes no further 
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constraints and we find the near horizon expansion of the fields to be 

F = /i(r-r+) + /2(r-r+)2 + ... 
D = di{r - r+) + d2{r - r+)^ + ... 
H{r) =Ho + Hi{r - r+) + H2{r - r+f + ... (4.1) 

'f{r) =<fo + <fi{r- r+) + (f2ir - r+f + ... 
p{r) = po +pi{r - r+) + p2{r - r+)^ + ... 

where r+ is the Schwarzschild radius of the black hole. By inserting these into the 
field equations and expanding order by order, appropriate boundary conditions can 
be found. This procedure leaves us apparently with four independent field variables 
at the horizon: /i, Hq, pq and ipQ for each choice of z or m. However, note that / 
can be shifted by a constant at the price of rescaling t, thus /i is essentially a gauge 
degree of freedom, which is tuned to achieve F — )■ 1 at infinity. Note also that the 
metric and field equations are invariant under the rescaling 

r ^ Ar , ^ ^ ^ , ^' ^ T' ^^-^^ 

which means that we are free to set the Schwarzschild radius of the black hole, r+, 
to 1. We choose, however, to keep r+ explicit in our expressions for clarity. 

Numerical solutions can be found by fixing either rh (for adS) or z (for Lifshitz) 
and using a shooting method to integrate out from the horizon, tuning the inital 
data to give a regular asymptotic solution. 

4.1 AdS Black Holes 

As a warm up, consider first the asymptotically adS solutions. For numerical sim- 
plicity we focus on the case where there is only one unphysical growing mode at 
infinity, i.e. rh G [a/6/5, v^lO — 2]. Note that in this range the exponents of 99 and H 
are complex. Via this process we find a two parameter family of asymptotically adS 
black hole solutions for a fixed value of rh. In light of the previous analytic findings 
these correspond to some combination of gauge and scalar charge. A priori we are 
free to choose any two of Hq., ipo and po as our free parameters and we shall choose 
them to be ipo and po- Figures | to |^ show examples of the asymptotically adS black 
hole solutions with m = 1.105. 

Figure § shows three solutions where p = 0, in which the black holes have only 
scalar charge. One has only a small scalar at the horizon, and the other two a more 
substantial dilaton charge, one positive and one negative. The small dilaton pertur- 
bation leads to a smaller perturbation in H and largely leaves F and D unchanged, 
which is consistent with the approximations made in finding the analytic solution 
( p.l7| ). Turning ipo up to 1 shows how positive dilaton charge reacts on the geometry. 
The perturbation in H grows and F and D are no longer equal, however, all three 
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Figure 5: Field profiles for asymptotically adS black holes with m = 1.105 and no 
i?— charge. In the left figure, the black lines correspond to the F function and the grey to 
D. On the right plot, the purple lines correspond to the dilaton, and the pink to the H 
field. In both plots the solid lines are with a small perturbation, (pQ = 0.2, the dashed and 
dotted plots to a larger dilaton charge, with ifo = l and ipo = —0.52 respectively. 

functions remain nionotonically increasing outside the horizon. Something more in- 
teresting happens however if we try to lower the value of the scalar at the horizon, 
i.e. letting (po < 0. Overall, very little happens to the geometry, however, as can be 
seen from the plot for (po = —0.52, the D and H functions cease to be monotonic 
and all the fields relax to their asymptotic values significantly more slowly than for 
positive scalar charge. Indeed, there is a critical value of (fo ~ —0.53, below which 
the charged black hole solution ceases to exist. This is because there is a runaway 
behaviour in V{ip) for cp too negative. This critical value shifts towards the origin as 
we turn on 5— charge, as can be seen by looking at the source for the if equation of 
motion at the horizon. 

Figure ^ shows solutions with 5— charge, while keeping ipQ = 0. Again, we show 
the comparison between a small and larger gauge charge. The small perturbation in 
P leaves the other fields largely unchanged, and comparing the analytic approximate 
and numerically generated profiles of P for the same value of m, one finds that the 
two appear identical to the naked eye. Increasing pq alters all the other fields as seen 
in the dashed plots, and once again F and D increasingly differ. 

The effect of both scalar and gauge charge on the black hole is shown in figure 
1^. Here, the dotted line shows a nearly critical negative scalar charge black hole, in 
which the 5— charge has been turned to near extremality (the temperature of this 
black hole is 0.024r+). 

4.2 Lifshitz Black Holes 

We now turn to black hole solutions that asymptote the Lifshitz spacetime defined 
by ( |2.11|) . A crucial difference between this case and the asymptotically adS case is 
that the background 2-form gauge field is now nonzero, p ^ 0- As a result we were 



-18- 







5 
log r 



10 



5 
log r 



10 



Figure 6: Field profiles for asymptotically adS black holes with rh = 1.105, with (pQ = 0. 
As before, the black, grey, pink and purple lines correspond to the F, D, H and ip fields, 
with the p-field being plotted in brown. The solid lines are for a small perturbation, 
Pq = 0.1, and the dashed to a larger charge, po = 1- 




Figure 7: Field profiles for generic adS black holes with rh = 1.105. The solid lines 
correspond to horizon data ipQ = 4\/2/5, pq = 1; the dashed lines to (pQ = — \/2/5, po = 1; 
and the dotted lines to (/?o = —0.2, po = 1-5. The colour coding of the plot is the same as 
figure ^. 

unable to find an exact analytic expression for the black hole similar to ( |3.9D , only 
perturbative or asymptotic solutions. We must therefore rely solely on numerical 
results for full solutions to the equations of motion. Another difference with the 
adS case is that all of the fields necessarily participate in the black hole solution. 
That this is the case can be seen by checking the eigenvectors of the perturbations 
around the critical point. All of the eigenvectors have (different) combinations of 
scalar, gauge and geometry components. Thus we expect that, for Lifshitz black 
holes, nontrivial profiles of all of the fields will generally be present. 

The numerical solutions are found in precisely the same way as in the adS case 
with the only difference being that the parameters of the theory are now defined 
by (|2.11| ) where the dynamical exponent 2; > 1 is used to fix the theory as opposed 
to m. For simplicity we present solutions corresponding to the upper sign choice in 
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( p.ll ) as this joins the adS branch and is more stable to integrate. By integrating 
the equations of motion ( |3.2|) to (|3.6|) we once again find a two parameter family of 
asymptotically Lifshitz black hole solutions for each value of z to which we assign the 
free parameters yjo and po- Using intuition from the adS case we suggest that these 
parameters relate to the scalar and 5— charge of the black hole, however, since all 
the fields participate in any asymptotic fall-off to the Lifshitz spacetime, this relation 
will not be completely straightfoward. 
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Figure 8: A selection of plots showing the field profiles for asymptotically Lifshitz black 
holes for z = 2. The colours are the same as for ads, with F, Z), H^ ip and "p = p/b being 
displayed in black, grey, pink, purple and brown respectively. The upper pair of plots 
explore the effect of varying c/?o, and the lower plots the impact of changing po. In each 
case, the labelling of the curves is defined by the initial conditions in the right hand plot. 

In all of the Lifshitz plots, we renormalize the z-dependence of the gauge field 
by plotting p = p/b, so that unnecessary variation with z is scaled out. Figure H 
explores the impact of varying the gauge and scalar initial conditions on a, z = 2 
Lifshitz black hole. The plots are reasonably self-explanatory, exploring the impact 
of altering ip^ (upper) and po (lower) relative to the fiducial black hole solution 
shown in each case by solid lines, whose horizon values of the dilaton and 5— field 
are the same as the asymptotic values. Most of these black hole solutions (see 
also figure |^) have extremely strongly warped geometries near the horizon, with the 
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"Newtonian potential", F, rising very sharply to a rather high maximum before 
falling to its asymptotic value of 1 from above. For example, in the fiducial solution, 
the maximum of F is around 30 (in units of r+), and lowering the horizon value of 
Po only exacerbates this effect. 

In the upper plots, showing some sample solutions for different (/jq, we see that de- 
creasing ^pQ below zero damps the F potential, which now looks more like a canonical 
black hole, with the example shown in the dotted plot, y^o = —0.5 being represen- 
tative of a nearly 'extremal' black hole, in the sense that the solution will cease to 
exist if if is lowered further (for the same reason as in the adS case) and also in 
the sense of the temperature dropping to zero. Increasing lpq on the other hand has 
the opposite effect, with 999 = 0.5 increases the sharp peak of F, which now has 
a maximum of around 80 (the dashed lines). Interestingly however, this variation 
with the dilaton horizon value is not monotonic, and as ^po is increased further, the 
amount of warping peaks, then subsides, and as we will see in the next section when 
we consider thermodynamics, this behaviour is mirrored in the temperature of the 
black hole dropping to zero. 

In the lower pair of plots, which explore changing the horizon value of the gauge 
field, we see that increasing po rapidly restores the F field to a more canonical form, in 
the case of the dashed (po = 1-5) and dotted (po = 2) plots. Since the temperature of 
these black holes drops, increasing the horizon value of p can be seen to be analogous 
to charging up a Reissner-Nordstrom black hole. Correspondingly, as can be seen in 
the dot-dashed plot, dropping po below its asymptotic value causes the black hole to 
become more strongly warped, and as we will see, hotter. 

Altering both ip^ and po produces a combination of these effects: increasing ipo 
first moves the black hole away from, then towards, 'extremality', and decreasing po 
always moves the black hole away from extremality. Correspondingly, the maximal 
value of 5— charge (tracked by po) first increases, then decreases, as we increase the 
dilaton charge ipo. 

Finally, the plots in figure ^ show the effect of altering the dynamical exponent, 
z, for a system with ipQ = 0.5 and po = 1.5. The plots show the field profiles for 
z = 2, z = 3 and z = 5.75, where the last value was chosen to be within the range for 
which all the eigenvalues of figure ^ are real. These plots show how z can alter the 
field profiles, particularly F and D, and that increasing z hastens the convergence 
to the Lifshitz solution. This was to be expected since, as can be seen in figure ^ 
increasing z largely reduces the eigenvalues governing each fields approach to the 
Lifshitz fixed point. 

4.3 Uplifting to Type IIA in 10 dimensions 

With the 6D solutions in hand, it is straightforward to uplift them to configurations 
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Figure 9: Plots showing the field profiles for asymptotically Lifshitz black holes with 
(fo = 0.5, po = 1.5 at different z. The black, grey, pink, blue and red lines correspond to 
F, D, H, p and (p respectively. The solid lines correspond to z = 2, the dashed lines to 
z = 3, and the dotted lines to z = 5.75. 



in Type IIA massive supergravity. Following O, P], we define 



V3m/ 

C(p) = cosp, S(p) = sinp 

A{p) = X C^ + X-^ S^ 

U{p) = X-^ S^ - 3X2 ^2 ^ 4^-2 ^2 _ g^-2 



(4.3) 
(4.4) 
(4.5) 
(4.6) 



as well as the constant k = (3ni5f^) /2. We can then write the ten dimensional, 
uplifted configurations as: 



dsl, = S'/''X'/' 



A'/'iLtBH, X n^) - 2k-^A^/^X^ dp" - \-'' A-^'I^X'^^ Y.'^h^'^f 



V2 ,_. 



F4 = — k-^ S^/3 C=^ A-2 UdpAe3 + V2 fc-^ S^/=^C X^ ^g ^3 A dp 
6 

_^ k-' S^/3 c fP ^ /,(3) A t/p + ^ k-^ S^/3 ^2^-1 ^-3 ^(3) ^ ^(1) ^ ^(2) 

V2 4^2 

+72 A;"^ S^^^C^A-^X-^dX A eg , 
A;2 



(4.7) 



G3 = 2^2— S'/^ G3 , F2 

g$^g-5/6^1/4^-5/4^ 

where 



0, 



i« 



with cr^*-* the left-invariant 1-forms on S^, and es = h^^^ A /i'^^-' A /i'-^-'. The parameters 
of the 6D theory are related to the Type IIA mass parameter via m = (2 m (7^/27) 
The uplift gives us some insight into the kinds of sources in lOD that give rise to 
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the adS and Lifshitz black holes in 4D. Notice that the ten dimensional RR F2 field 
strength vanishes, while the RR F4 field and the NS G3 field are switched on in 
several directions. The main difference between the black hole solutions discussed 
here and their pure Lifshitz/adS counterparts presented in [^] is in the non-trivial 
r-profiles, which imply that F4 contains an additional component in the direction 
dX A £3. The configurations could be interpreted as a system of D-branes and NS- 
branes, analogously to [Q. 

5. Thermodynamics 

Having determined the black hole solutions for our supergravity set-up, it is interest- 
ing to investigate some of their general properties besides their field profiles. In this 
section we make the first few steps of this exploration by studying the dependence 
of the temperature of these black holes on the initial parameters ipo and po, as well 
as commenting on entropy. 

The temperature of the black hole is given by 



-1 

Air 



T= ^VD'F' \r=r+, (5.1) 



which is calculated directly from the numerical solutions. Figures ^ and show 
the temperatures of both asymptotically adS, and asymptotically Lifshitz black holes 
as a function of our initial parameter pQ and ipQ. In all the temperature plots, we 
have shown the temperature normalized at r+ = 1 for simplicity. We discuss the r_|_ 
dependence of the temperature at the end of this section. 

In figure 0, plots are shown of the temperature of an asymptotically adS black 
hole with rh = 1.105. On the left, the plot is shown as a function of the 5— charge 
for a range of ipo, on on the right as a function of (po for a range of Pq. The left 
plot shows the expected behaviour of a charged black hole, in that adding charge 
reduces temperature monotonically to zero at an extremal limit. The effect of scalar 
charge in this case is more interesting. At zero -B— charge, the impact of increasing 
(fo is to increase the temperature of the black hole, and one might expect therefore 
that the allowed 5— charge range is increased. However, an interesting phenomenon 
occurs. As the black hole becomes more and more charged under the scalar field, the 
maximal amount of 5— charge we are able to add starts to drop, and at very high 
scalar charges we can no longer add much gauge charge. That the "extremal" limit 
should not be a simple sum of the two charges, but some more complex combination 
is an interesting difference from most black holes with more than one charge. 



Figure |TT] shows the corresponding plots for a Lifshitz black hole, with z = 2 
taken as an example. As with the adS black hole, increasing po reduces the tempera- 
ture, and once again, we see that there is a finite range of ipo for which the black holes 
exist. However, because the Lifshitz spacetime has a nonzero background i?— field. 
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Figure 10: Plots of the temperature of an adS black hole with m = 1.105 shown as a 
function of po (left) and (po (right). In the left plot, the dashed pink data in inverted 
triangles corresponds to (po = —0.2; the black data in dots to fo = 0; the dashed blue with 
squares to (/?o = 1; the dotted red with diamonds to c/?o = 2; the grey triangles to 939 = 3; 
and the dashed black with open circles to (po = 4. In the right plot, the lines run from 
Po = in black with dots at the top, to po = 4 in dashed purple with inverted triangles at 
the bottom in increments ApQ = 1. 

there is a clear difference in the temperature as a function of po — ^ 0. For the Lifshitz 
black hole, the temperature increases sharply as we reduce the initial value oi po, and 
would appear to diverge as p — )■ 0. In the absence of analytic arguments we cannot 
say definitively that T diverges, however, our numerical integrations become more 
and more extreme as we reduce pQ. It is worth noting that we have presented our 
temperature plots renormalized to r_|_ = 1, clearly, dropping r+ drops the tempera- 
ture, so it is possible that we can achieve po ~^ by taking r+ — )■ 0, indeed, such a 
spacetime would represent a flow from a Lifshitz space in the UV to an adS space in 

m 



the IR, 



However, the results of ||2^ would indicate that this will only happen 
for a specific value of ipo, namely, the one corresponding to an adS solution. Thus, 
we would expect the generic po ~^ limit to be singular. 

In figure we see how altering z affects the temperature of the black hole by 
looking at the variation of temperature with po for sample values of z, and exploring 
in detail the z-dependence for sets of representative initial data. In general, we 
see that increasing z raises the temperature for small po but lowers it for larger pQ, 
although for temperatures close to 1, the temperature seems to first increase then 
decrease with pQ. 

Finally, we consider the entropy of the black holes, which can be computed from 
the area of the horizon as: 



S = -rlH{T, 



(5.2) 



per unit volume (and setting Gg = 1). Figure ^ shows how the entropy of a 2; = 2 
Lifshitz black hole varies as a function of temperature at fixed r+ = 1 for different 
values of initial data. The two plots show set of data for S{T) having fixed (^0 and 
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Figure 11: Plot of the temperature of asymptotically Lifshitz black holes for varying po 
and ifo, with z = 2. On the left the plot depicts the temperatures of z = 2 black holes as a 
function of po for ipo = in black with circular data points, (po = 1.5 in dashed blue with 
square data points, and for (/Jq = 3 in dotted red with diamond data points. On the right, 
the temperature is shown as a function of ipo for po = 1 in black (circles), po = 1.5 in blue 
(dashed/squares), and po = 2 in red (dotted/diamonds). 
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Figure 12: Plot of the temperature of asymptotically Lifshitz black holes for varying poi 
ifo, and z. On the left, the curves give the temperature of a Lifshitz black hole with (pQ = 1 
as a function of po for different values of z: black (circles) to z = 2, blue (dashed/squares) 
to z = 3, and red (dotted/diamonds) to z = 5. The three curves appear to intersect at a 
single point, however, the resolution of the data is insufficient to be sure if this is exact. 
On the right, the temperature is plotted as a function of z for initial data {^o,Po) = (0, 1) 
in black (circles), {^o,Po) = (1;2) in blue (dashed/squares), and (fOiPo) = (0, 1.25) in red 
(dotted/diamonds) . 



varying po (left) or data for fixed po and varying (po (right). In the left plot, we see 
that just as varying po at fixed cpo has a much more uniform effect on temperature, 
so varying po at fixed ipo has a somewhat more consistent effect on entropy, although 
curiously the entropy generally drops as we increase the temperature. For ipQ = 
however, the entropy remains fairly constant. Since we have fixed r_|_ = 1, we would 
not necessarily expect the entropy to vary hugely with po, at least by analogy with 
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the Reissner-Nordstrom solution. 

When exploring the S{T) plot for varying ipo however, the picture becomes much 
more interesting. We would expect entropy to vary much more strongly with cpo at 
fixed r+, since we have seen from our eigenvalue analysis of the perturbations that the 
two "scalar" modes, the dilaton and the internal breather H, are very much coupled 
by the equations of motion. We therefore expect that altering (pQ will alter H{r^) 
and hence the entropy to a much greater extent, and this is indeed what we see. 
However, what is interesting is the modulating behaviour of both the temperature 
and entropy as a function of (po. We see that at a given temperature and po, there 
are two possible values for the scalar charge, one with higher entropy that the other. 
Although it is not entirely clear from the plot, the curves have (po increasing in 
a clockwise direction, hence it is the black hole with lower ipQ that is entropically 
preferred. This indicates that these black holes will likely have scalar instabilities, 
perhaps shedding scalar charge to increase their overall horizon area. How this is 
consistent with the usual concept of a black hole accreting to increase its area might 
prove an interesting investigation. 

Finally, we should comment on the impact of varying r+ : Because of the scaling 
symmetry present in the equations of motion, all of the numerically computed fields 
are dependent on r/r^, and thus simply stretch with r_|_. In particular, the horizon 
value of H does not change with r+, and the derivatives of F and D just scale as l/r_|_. 
Thus, the entropy scales as r^, independent of the value of z, and the temperature as 
r^. (Of course, the entropy and temperature vary with the initial data of the charges 
as we have seen.) The variation of entropy with temperature is therefore explicitly 
the expected relation S oc T^/^ for a field theory in fiat 2+1 dimensions. 
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Figure 13: Plots of the entropy of a z = 2 Lifshitz black hole as a function of the 
temperature with r+ = 1. On the left, the data points represent different values of po 
for ifQ = 0, 1,2, and 3, in black (dots), dashed blue (squares), dotted red (diamonds) and 
dot /dash purple (triangles) respectively. On the right, the data points represent different 
values of cpo for po = 0.75, 1, 1.5 and 2, in black (dots), dashed blue (squares), dotted red 
(diamonds) and dot/dash purple (triangles) respectively. 
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6. Conclusions 

In this paper we have buih string black hole configurations that asymptote Lifshitz 
geometries for general dynamical exponents z > 1. We began by finding these solu- 
tions in a supergravity theory that corresponds to a consistent truncation of massive 
type IIA string theory, namely the maximal gauged six dimensional supergravity of 
Romans, and then uplifted them to massive type IIA supergravity. They are related 
to asymptotically adS black hole solutions in the same supergravity, also studied 
here, and extensions of the pure Lifshitz geometries found in 0. 

Our supergravity theory has rich field content, and yet allows simple pure adS4 x 
H2 and Li4 x H2 compactifications. By expressing the full set of equations in the 
form of an autonomous dynamical system, we studied perturbations around these 
exact solutions, determining analytically all the possible asymptotic behaviours of 
static, planar adS and Lifshitz black holes for our theory. We found that various 
asymptotics are allowed, depending on which fields are switched on and on the values 
of the parameters involved, in particular the dynamical exponent z. 

Starting with the adS case, one exact analytic solution is of course known: adS- 
Schwarzschild. By studying this background, and exciting the supergravity two-form 
potential B2 and dilaton in the probe limit, we acquired insight into the properties 
of adS black holes. In particular, apart from the horizon size, the black holes are 
characterized by two quantities, ipo and po, which we are motivated to call scalar and 
S-charge. After exploring the adS case, we applied the intuition gained from it to 
analyze Lifshitz configurations, for arbitrary dynamical exponent z > 1. In fact, a 
non-trivial scalar and 5-charge is a necessary ingredient for obtaining asymptotically 
Lifshitz geometries. Although the notion of a charged black hole in the Lifshitz 
case is less straightforward than for adS, since any Lifshitz configuration is always 
characterized by the presence of all fields, we thus determined the independent, 
tunable quantities characterizing Lifshitz black holes. 

We were able to formulate a precise analytical understanding as to how asymp- 
totically Lifshitz configurations reduce to adS ones, in the limit ;z — )■ 1, which could 
be helpful in the future for comparing predictions of field theory duals of adS and 
Lifshitz theories. Then, we studied numerically how the asymptotics match black 
hole horizons at finite values of the radial coordinates, computing two-parameter 
families of black holes, for both adS and Lifshitz asymptotics. 

One common feature displayed by the Lifshitz black holes is the presence of a 
non-monotonicity in the time and radial metric potentials, F and D, particularly well 
illustrated in figure p. This peaking of the Newtonian potential becomes extremely 
marked at lower ipo and po? with F peaking at over 50 or higher. Although it 
might appear as if a large peak in F could provide a barrier to approaching the 
event horizon, in fact this is outweighed by the strong Lifshitz r^^ warp factor, and 
particles are always pulled into the black hole. 
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As well as studying the field profiles for the black hole solutions, and how they 
vary with S-charge, scalar charge and dynamical exponent 2;, we initiated a study of 
the thermodynamics of our black holes. We find the 5-charge to play an analagous 
role to the charge in the Reissner-Nordstrom black hole, with temperature decreasing 
as the charge increases, however, the zero charge limit is more subtle for the Lifshitz 
black holes. Meanwhile, adS/Lifshitz black holes exist only for a finite range of 
scalar charge, due to the Liouville type scalar potential, and in general for a given 
temperature there are two possible scalar charges, perhaps surprisingly the smaller 
one being entropically preferred. The interplay between the scalar and i?-charge is 
also interesting, in particular a high scalar charge reduces the amount of possible 
S- charge. 

Since we are interested in holographic condensed matter applications, we have fo- 
cussed on planar black hole geometries, corresponding to a finite temperature bound- 
ary theory in fiat 2+1 dimensions. In simpler models, with pure geometry planar 
black holes, the scaling symmetry renders black holes with different horizon sizes, and 
thus temperatures, physically equivalent. However, in our supergravity construction 
the necessary, additional scales of 5-charge and scalar charge introduce the possibil- 
ity of phase transitions. It would be interesting to investigate this possibility further, 
as well as their interpretation in the dual field theory^. Moreover, the extra fields 
present in the supergravity theory, which we have chosen not to activate, represent 
other possibilities for interesting black hole charges and phase transitions. 

Another aspect that deserves further thought would be the brane interpretation 
of our type IIA configuration. In the end, we hope that having explicit string theory 
embeddings of Lifshitz black hole geometries will help to develop the holographic 
description of interesting anisotropic condensed matter systems. 
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'''Phase transitions in a phenomenological model that strongly resembles our stringy supergravity 
setup have been studied very recently in p5| . 
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A. The autonomous dynamical system 

In this appendix we reformulate the system of supergravity equations of motion as 
an autonomous dynamical system. This requires the use of a different gauge to 
the numerical work, however the translation between the gauges is straightforward. 
Choosing d = — 0/v2 — ln-L(yf, in the general planar metric ( p.l3| ), ensures that only 
negative exponents of (p appear in the equations of motion. Then, defining 



V ^ 9 

X5 = 2/i' + y20' , X6 = 2c' , X7 = e'2=+^^P , Xg = e'^^+^^^P' (A.l) 

yields the dynamical system 

X[ = -X1X3 X', = -X2X5 (A.2) 

^3 = -x',+ x^ (xi-^]-U xl - ^-xt 



A J A 'A 
-Xs[^^ + X, + X,yf (A.3) 

Xf\ , 1 , ^2 xf 



X', = X^,+X^(3Xl + ^y\ + Xl 



-^4(^^ + X, + Xe)+f (A.4) 
^5 = (^?-l)(2X,^ + ^)-^ + ^ + ^4X6 + ^ 

HX. - 2X3) (^ + Xe) - ^l + ^f^ + ^l (A.5) 

X'e = X!-X^^Xl + '-^^+\ + X^,-^ 

-Xa(^i^ + X3 + Xe)-f (A.6) 

X; = X8-X6X7 + X3X7 (A.7) 

X^ = -X8(ix4 + X5 + ^) + X7 [xf + 4X2^] (A.8) 

in which the solution lies in an invariant 7D submanifold described by the constraint 



2X4(X6 + X5 - X3) + 4(X5 - X3)X6 + Xi + Xi - 2X5X3 - Xi - Xi 
= 4:^X2 - 4X2 _ x2(4X2' + Xf) + {1 + AXl - X\) . (A.9) 

For a critical point, ( |A.2| ) implies X3 = X5 = 0, and ( |A.7| , |A.8| ) give: 

Xs = XgXr (A.IO) 

X4X8 = 2X7 [X^ + 4X2] (^j^j-^^ 
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These are solved by either Xj = Xg = 0, or 

X4X6 = 2X^ + SXl . (A. 12) 

We can distinguish two cases: 
• Case 1: adS 

In this case X7 = Xg = 0. Solving (|A.3| - |A.6|) gives 



Xl = l-\ll-\xl (A.13) 



with the constraint 



^2-g^6 4 + 2V 

x, = x. 


l.-\n 


'^:-sii-xi>-2, 


/-!-« 



(A.14) 
(A.15) 

(A.16) 



selecting two values for Xg for each charge. In terms of the supergravity parameters, 
g = 2/Xq, rh = 2x\/Xq. To get the solution in the original area gauge, note that 
r = e^ = e^'^^l'^ . These critical points form a curve of adS solutions in the phase 
space, with the curve intersecting the invariant submanifold at two points in general. 
To analyze the nature of the critical points, one takes small perturbations, and 
finds the eigenvalues and eigenvectors of the perturbation operator matrix given by 
(5X^ = Mo/^X^, i.e. Mt>('^ = Ajf*^*\ There will always be one zero eigenvalue to 
this matrix, corresponding to moving along the adS solution curve. The remaining 
7th order polynomial can be factorized, yielding the eigenvalues plotted in figure 0, 
although note that to obtain figure 0, we have transformed our coordinates to the area 
gauge to get the Aj to correspond to the r fall-off exponents: Aj = 2Aj/X6 = g\i. 
• Case 2: Lifshitz 
Here Xg = XgXy, and ( |A.12| ) holds, and solving the remaining equations gives: 




X4 = X6(l + X^) = J ' , ' (A.17) 



2 (5 + X|) T v/2(5 + X|) 
^' - (1 + X|)(5 + X|) ^^-^^^ 

^, _ (X| + 3){X^ - 2) ± 2v/2(5 + X|) 

4(5 + X|)(l + X|)2 • ^'^■''> 



For X^ < ^2 - (1 + a/17 + 4a/2)/2 ~ 3.3, in ( |70^ ) only the upper branch choice 



gives a real solution for X2. Finally, the constraint determines the charge: 



(X| + 5) ((X| + 3)(X2 - 2) ± 2v/2(X2 + 5) 

qg = -^ , . (A.20) 

3(X2 + 3)t2v/2(X2 + 5) 
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Clearly X| = z — 1, and this is equivalent to the exact solution ( p.ll| ) in the area 
gauge. 

The analysis of the perturbations around the Lifshitz critical points, although 
conceptually identical to adS, is algebraically more involved. The eigenvalues pair 
around —{z + 2)/2 (in area gauge) with the variance given by the square root of 
the solution of a cubic equation. While this can be written in closed form, it is a 
rather long and unilluminating expression. Figure § shows a plot of the eigenvalues 
renormalized for a r fall-off: Aj = a/2z(4 + z)Aj. All of the eigenvectors of the 
perturbation operator have nonzero components in the vector, scalar and geometry 
directions. 



B. Exact Lifshitz solutions 

So far in the literature it has been possible to obtain exact Lifshitz black hole (LiBH) 
solutions only in some phenomenological models, where the matter content is engi- 
neered to support the desired geometry. In this appendix we give a brief account of 
these approaches, and extend them to dilatonic models, which may be more easily 
embedded into supergravity and string theory. We then show that such simple ana- 
lytical solutions cannot be found in the Romans' 6D supergravity that is the main 
subject of this paper. 

Analytic LiBH solutions have been found for essentially two types of 4D Einstein 
gravity systems (see also |]2D| for other possible extensions). The first [AAAm) 
contains, besides gravity, a cosmological constant, a massless abelian gauge field J-2 



and massive abelian gauge field F2 with mass m |T^. In 4D, the massive gauge field 
is equivalent to a 2-form B2 and a massless gauge field F2 with non trivial Chern- 
Simons terms F2AB2, as studied in ^ . The second system ( A^0) consists of gravity, 
a cosmological constant, a number of massless abelian gauge fields J-'g and a massless 
scalar field with dilaton-like couplings to the gauge sector |^, 18|. An obstruction 



to straightfowardly embedding these setups into supergravity and string theory is 
their absence of a genuine dilaton field. For example, many supergravity theories, 
like Romans' 6D supergravity, necessarily contain a dilaton field in the supergravity 
multiplet. Indeed, it seems unlikely that the AA(f) system could be embedded in 
string theory, without generalizing the cosmological constant to a genuine dilatonic 
potential. It is then easy to check that in the presence of a non-trivial dilaton 
potential, the dilaton equation of motion prevents the Lifshitz asymptotics. This 



can be observed from the general solutions discussed in ^7_ . 

The AAAm system is more interesting. Building on [^, |], a similar model, 
though without the massless vector and with two additional dilatonic/radion scalars. 



was obtained in |21| via a consistent massive truncation of Type IIB on an arbitrary 



Einstein space times S^, and used to derive numerical stringy LiBHs. Here, we gen- 
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eralise the system to a general dilatonic theory, with generic dilaton-matter couphngs 
and dilaton potential, and search for analytical LiBHs. 

B.l 4D LiBHs with constant dilaton 

We consider the four dimensional case, which is sufficient to illustrate our strategy. 
As well as massless and massive abelian gauge fields, we add a dilaton field with 
couplings A, 0" to the gauge fields and a general potential V{(f)). In this section 



we follow closely the discussions in [|1^, |I^, and assume their mostly plus metric 
signature conventions. 

We take the general action: 



S 



d^Xy/g 



1 p-^0 






(B.l) 



with corresponding equations of motion: 



(B.2) 
(B.3) 

--d^ [V^g^'dA] = ^ - —^F^ 7. A —^ (B.4) 



-9 



V 



-At/) 



Ri^v = ^d^cpd^cp + —g^^ + ^ 



2 

2F'F,,-^F' + 2m'A,A 



fj,-' '-V 



2J-^ J-^\ —J- 



Consider now the metric Ansatz 



ds' 



r''h{r) dr + 



111 
— — — — - + r aXj , 



(B.5) 



(B.6) 



sourced by J>i, FxiX2 ^^^d F^t, plus a constant dilaton field (j) = const. Solving the 
field equations for the forms gives: 



•Frt — Qi^ 

dr[r^-'drAt] 



■J X\X2 ^1 



m^At 



(B.7) 
(B.8) 



From the Einstein equations R\ — R'^, we obtain the solution for the gauge field: 



and then using 



m 
we find the metric function: 

^^ Cl 1 ^2 



ni 
h = — + 



2z {z-2)r- 



+ 






(B.9) 



(B.IO) 
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for 2 7^ 2, and an analogous expression for z = 2. 

Now the dilaton and remaining Einstein equations can be solved provided that 
the following constraints among the parameters are satisfied: 



z = A 

-—[2{z + 2) + z{z-l)] = - 
2X{z - 1) C2 = ^e^'^HQl - Ql) 



Co<0 



X{z-l)m'[^ + -^ 



dV 



(B.ll) 
(B.12) 

(B.13) 

(B.14) 
(B.15) 

(B.16) 



The first four constraints above correspond to those in ^Tj where there is no dilaton, 
and the dilaton adds two more. In the end there are four non-trivial constraints on 
the four solution parameters C2, Qi, Q2 and (p, plus the mass parameter m and any 
gauge couplings that appear in V^cj)). So, provided that these constraints can be 
solved consistently, we can avoid tuning the dilaton couplings A, a. 

B.2 6D Romans' LiBH with constant dilaton? 

Romans' 6D supergravity has a strong resemblance to the dilatonic theory we just 
discussed, or its Chern-Simons equivalent. Besides having two extra dimensions, the 
main difference is that the 2-form potential has not only a Chern-Simons term but 
also a mass term. Nevertheless, it is straightforward to apply the above strategy 
to search for LiBH solutions with constant dilaton in Romans' supergravity. These 
solutions would be orthogonal to the ones we discuss in the main text, since they 
involve turning on additional fields. 

Indeed, whereas in the main text we activated only one of the gauge fields, we 
now add a non-trivial configuration for the second gauge field, taking the Ansatz 
(PD plus: 

J'rt = Qir'~' T.,., = Q2. (B.17) 

The expression for the metric function that follows is: 



h{r) 



m.Q2e' 



-2^200 



2)v^ 



Ir^ 



r^ 2z 



^2g-3V20O+4^e-V2<^O 



(R18) 



for z y^ 2, and an analogous expression for z = 2. Unfortunately, the constraints 
coming from the field equations can only be solved by a pure Lifshitz configuration, 
with Qi = = Q2 and Ci = 0, and the relations ( p. lip . Thus we cannot construct 
a simple analytic solution via this method. 
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